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1. Background 

Through recent observational programs, such as 
FIRE, the structure of the cirrus cloud layer has re- 
cently been determined to be quite complex. Work 
of Sassen, et al. (1989) with lidar, aircraft and radar 
found a wide range of spatial scales from generating 
cells of about 1 km in size to mesoscale complexes of 
these cells on the order of 10’s km, which are imbed- 
ded in the larger cirrus cloud. Smith, et al. (1989) 
with aircraft and Grund and Eloranta (1988) with li- 
dar also found much spatial variability in the cirrus 
cloud parameters. 

This research in studying the spatial scales of the 
cirrus, used data collected during the flight legs of 
the NCAR Sabreliner aircraft on four days during 
the FIRE Cirrus IFO to study the spatial scales of 
the cirrus, and will concentrate on the scales of the 
horizontal wind. The days examined are 22, 24, 28, 
and 31 October 1986. A variety of conditions ex- 
isted on those days. On the 22nd and the 24th, a 
cirrostratus deck with convective cells was sampled 
in the vicinity of the exit region ridge of a closed low 
pressure system. Starr and Wylie (1988) classified 
the cirrus as “closed low cirrus”. The cirrus of the 
28th was classified as “warm front cirrus” (Starr and 
Wylie, 1988) because the cirrus was found ahead of 
the upper level short wave trough, capping the short 
wave ridge. The cirrus was optically thin and ten- 
uous. The cirrus of the 31st was classified as “cold 
front cirrus” (Starr and Wylie, 1988) or “jet stream 
cirrus” (Conover, 1960) and was found on the eastern 
side of a short wave trough. A jet with wind speeds 
of about 50 m/s was measured in the cirrus top. 

2. Data Analysis 

The spatial scales of the cloud features can be de- 
scribed by power spectra (or spectral density graphs) 
and cumulative variance graphs. The cumulative 
variance graphs were created by first using a Fast 


Fourier Transform (FFT) to create variance spectra. 
The variances were then summed in a cumulative 
fashion from the largest scalelengths (wavelengths) 
to the smallest. No detrending was done to the orig- 
inal data, and no smoothing or averaging was done 
to the spectral points. All the spectral points were 
included. This means that the values of the first five 
to ten spectral points of the large scalelengths should 
only be considered to be qualitatively correct. The 
cumulative variance at smaller scalelengths should be 
correct because a more accurate representation of the 
variance at the larger scalelengths should only redis- 
tribute the energy amongst the larger scalelengths. 
Note that each scalelength (or wavelength) contains 
both a “peak” and a “valley”. The following graphs 
of cumulative variance may be interpreted as showing 
the percentage of variance explained by scalelengths 
larger than the given scalelength. 

The power spectra likewise were created using an 
FFT and differ from the variance spectra by only a 
constant factor. However, the data for the power 
spectra were detrended and the first five spectral 
points were not used. The power spectra were 
smoothed with a six point average. 

The optical properties of a cirrus cloud resulting 
from the downwelling shortwave and infrared radia- 
tion can be expressed by the infrared downward effec- 
tive emittance (c* j) and the shortwave effective ex- 
tinction U), which are defined by the following equa- 
tions. 
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where T is the mean cloud temperature, H T j is the 
mean cloud top incident shortwave irradiance, and 

J 

H t | is the mean cloud top downwelling infrared 


311 


PRECEDING PAGE BLANK NOT FILMED 


FIRE Cirru* iFO 25 OCT 56 



FIRE Cirrui IFO 28 OCT 86 



ScalVttnqtft (m) 


Figure 1: Cumulative variance graph of infrared Figure 2: Cumulative variance graph of shortwave ef- 

downward effective emittance. fective extinction. 


irradiance. This makes €* J, and £ dependent only on 
the down welling irradiance values {H§ w l t H^ R 1) 
at cloud base. Note that there may be a significant 
reflected component of upwelling irradiance inherent 
in e* 1 and 

S. Results 

The spatial scales of infrared downward effective 
emittance and the shortwave effective extinction are 
described by the cumulative variance graphs of Fig- 
ures 1 and 2 for the 28th. Most of the variance 
is found at the large scalelengths. For the infrared 
downward effective emittance, 90% of the variance 
is explained by scalelengths larger than 10 km; and 
for shortwave effective extinction, 90% of the vari- 
ance is explained by scalelengths larger than 5 km. 
Thus, most of the variability of these clouds, as seen 
by the Sabreliner radiometric instrumentation, is de- 
termined by the larger features. Also, the similar- 
ity between the shortwave effective extinction curves 
shows that the clouds have the same spatial scales 
even though Cloud Sample 2 is less optically thick 
than Cloud Sample 1. A higher percentage of the 
variance is found at the larger scalelengths in the 
infrared downward effective emittance than in the 
shortwave effective extinction. This is because the 
shortwave radiation, being highly anisotropic, is not 
smoothed as much by the hemispheric view of the 
radiometric instrumentation as is the infrared. The 
scales are then better represented in the shortwave 
effective extinction graph. The cumulative variance 
graphs of the 28th show this dominance of the larger 
scales in providing the majority of the variance, and 
are representative of the other days. 

Representative power spectra of the alongwind and 
crosswind components of the horizontal wind are seen 
in Figures 3 to 6, and of the corresponding down- 
welling shortwave and infrared irradiances in Figures 


7 to 10. The power spectra show a range of slopes 
from k 6 / 3 to k 3 . The wind power spectrum of the 
22nd shows a k 6 / 3 slope, while the wind spectrum 
of the 28th shows a A;” 3 slope. The power spectra 
of the wind for the 24th and the 31st have a range 
of slopes in between the Ar 6 / 3 and Ar 3 slopes. The 
power spectra of the irradiances also have slopes con- 
strained by the k 6 / 3 slope and the A: -3 slope, though 
the slopes are not always the same as the correspond- 
ing wind spectra. It is probable that the downwelling 
irradiances, modulated by the cloud elements, have 
their scales determined by eddies of the wind. 

4. Discussion 

The slope of the power spectra of the horizontal 
wind can be used to determine the energy transfer 
between the scales. Three dimensional turbulence 
has a k 6 / 3 slope in the power spectra showing an 
energy cascade to smaller scales or larger wavenum- 
ber (A;). The power spectra of effectively two di- 
mensional turbulence displays a AT 3 slope and a 
k 6 / 3 slope in different wavenumber regions. The 
k 3 slope is a result of the transfer of eddy enstro- 
phy (mean-squared vorticity) or the transfer of the 
square of quasi-geostrophic potential vorticity toward 
smaller scales. The Ar 6 / 3 slope is a result of an en- 
ergy cascade to larger scales. These cascades have 
been described by Gifford (1988), Charney (1971) 
and Kraichnan (1967). Gifford (1988) suggests that 
in the atmosphere the net energy cascade of the two 
dimensional AT 5 / 3 slope is to smaller scales. An alter- 
nate theory offered by VanZandt (1982) suggests the 
slope can be explained by internal gravity wave dy- 
namics. However, Gage and Nastrom (1986) did not 
find the wave model to be completely supported by 
the observations. They suggested that observations 
of the background vertical velocity spectrum supports 
the wave model while the turbulence model is sup- 
ported by the horizontal wind spectrum. It needs to 
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be noted that the two dimensional turbulence also 
known as geostrophic turbulence is normally applied 
at the synoptic scale and mesoscale motions. In this 
case the scales are smaller, having measured small 
mesoscale and large microscale motions (7 mesoscale 
and a microscale as characterised by Orlanski, 1975). 

The scales are likely too small for Charney’s 
geostrophic turbulence, however, the close proxim- 
ity of the tropopause and the strong temperature in- 
version is likely constraining the eddies with a fc” 3 
slope into an effectively two dimensional flow with 
an enstrophy cascade. The stabilities of the layers 
as described in Table 1 further support this hypoth- 
esis. Table 1 contains the vertical temperature gra- 
dient for the layer, the wind shear for the layer, and 
the Richardson number of the layer. The Richardson 
number is defined as 



d$ 


(£) 2 


where 6 is potential temperature, z is height, u is 
the wind speed, and g is the gravitational constant. 
The layer measured of the 22nd with the slope 

is the most unstable of case, while that of the 28th 
with the Jk -3 slope is the most stable, supporting the 
hypothesis that strong stability restricts the eddies 
to an approximate two dimensional flow. 

To summarize, the majority of the variance of the 
cirrus cloud properties is found at the larger scales. 
The radiative properties of the cloud appear to have 
their scales determined by the surrounding eddies, 
which modulate the cloud elements. The amount of 
variance or energy contained at each scale appears 
to be determined by the stability of the layer, which 
determines whether there is a two dimensional or a 
three dimensional turbulent cascade. 

Table 1. Stability parameters for the cirrus layer. 


Day 
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22 

1.9 

3.3 

5.1 

24 

2.5 

3.7 

5.3 

28 

2.9 

2.1 

19.2 

31 

3.1 

4.0-4.6 

4. 3-5. 8 
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Figure 3: Power spectra of the alongwind and cross- 
wind components of the horizontal wind: 22 Oct 


1986 . 
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Figure 4: Power spectra of the alongwind and cross- 
wind components of the horizontal wind: 24 Oct 


1986 . 
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Figure 5: Power spectra of the alongwind and cross- 
wind components of the horizontal wind: 28 Oct 


1986. 
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